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Results of our autocorrelation measurement performed on Fujitsu AP1000 are reported. We analyze (i) typical 
autocorrelation time, (ii) optimal mixing ratio between overrelaxation and pseudo-heatbath and (iii) critical 
behavior of autocorrelation time around cross-over region with high statistic in wide range of for pure SU(3) 
lattice gauge theory on 8 4 , 16 4 and 32 4 lattices. For the mixing ratio K, small value (3-7) looks optimal in the 
confined region, and reduces the integrated autocorrelation time by a factor 2-4 compared to the pseudo-heatbath. 
On the other hand in the deconfined phase, correlation times are short, and overrelaxation does not seem to matter 
For a fixed value of K(=9 in this paper), the dynamical exponent of overrelaxation is consistent with 2 

Autocorrelation measurement of the topological charge on 32 3 x 64 lattice at (3 = 6.0 is also briefly mentioned. 
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1. Aims of our study 

Aims of our autocorrelation measurement are 
to investigate the followings. 

1) typical autocorrelation time r 

2) operator dependence of r 

3) (3 dependence of r 

4) how overrelaxation(OR) reduces r or not ? 

5) optimal mixing ratio between OR and pseudo- 
heatbath 

6) lattice size dependence of r 

7) origin of large error of r 

8) critical exponent of r 

9) r for topological object on large(32 3 x 64) lat- 
tice 

Points 1 to 5 were previously reported upon at 



LAT92 §. Points 6 - 8 are new, as well as some 
preliminary results about 9. 

2. Notations and Parameters 

Notations and parameters in this paper are as 
follows. 

1) OBSERVABLE 

We adopt the lxl wilson loop on 2 lattice 
blocked from original lattice as an observable of 
our autocorrelation measurement. 

2) MIXING PARAMETER K 

To check the efficiency of overrelaxation, we intro- 
duce the mixing parameter K which means that 
the ratio between the Brown- Woch microcanon- 
ical updating]^] and Cabibbo-Marinari pseudo 
heat bath updatingra is set to be K : 1. 
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3) AUTOCORRELATION FUNCTION 

The autocorrelation function p(t) of observable O 
(=1x1 wilson loop) is defined as: 

(t) = (m)-QA)m+t)-o B )) (1) 

y/((o(i)-o A )*){(o(i + t)-o B ) 2 ) 

where ( ) denotes average over i and Oa — (0(i)) 
and B = (0(i + t)). 

4) AUTOCORRELATION TIME 

The autocorrelation time is defined in this paper 



as: 



N 



Tint = P (0) + 2j2p(t) 



t=l 



N-t 

N 



(2) 



where N is determined so that Ti nt is maximized, 
but N < 10 % of the total sample and N < 3rj ra t. 
5) LATTICE SIZE and (3 

8 4 : 29 /? values for X=9, 12 /3 values for X=0 
and 4 /3 values for i^=0-20 

16 4 : 45 /3s for X=9 and 2 /3s for #=0-12 

32 4 : 4 /3s for #=9 



3. RESULTS 

In Fig.l and Fig.2 we show the mixing param- 
eter dependence of the autocorrelation time r on 
8 4 and 16 4 lattice respectively. In both cases, 2 
/3 values are presented, just below the crossover 
corresponding to confinement (a), and above it 
(b). Note the expanded scale and short autocor- 
relation times in (b). In the confined phase, a 
small value of K ( 3-7 ) seems to reduce the au- 
tocorrelation time by a factor 2-4 compared to 
pseudo-heatbath. Above deconfinement, no clear 
K-dependence can be seen. 

t is affected by large errors in these figures. To 
check the origin of this error, we prepare 40000 
sweeps on 16 4 lattice at /3=6.40 after the ther- 
malization. We divide 40000 sweeps into 4 bins 
of 10000 sweeps each. Next we calculate p(t) in 
each bin. The result is shown in Fig. 3. We can see 
large discrepancy between the results obtained 
from these 4 bins. Clearly 40000 sweeps are not 
sufficient to obtain precise values for the autocor- 
relation function and time in this rather extreme 
case, very close to the deconfinement transition. 

Next we show the results of t on different lat- 
tice sizes in Fig. 3. In this figure we plot r on 8 4 , 



16 4 and 32 4 lattices as a function of the inverse 
lattice spacing a -1 . In order to convert the lat- 
tice coupling 0(= 6/g 2 ) to the lattice spacing a, 
we use results of ref.4, obtained by a Monte Carlo 
Rcnormalization Group analysis. From Fig. 4 we 
conclude that the critical exponent z defined by 



(— ^-) — ^ XaL y 

T a L CbL 



(3) 



is consistent with 2. A detailed analysis is now in 
progress. 

Finally we mention briefly our measurement of 
topological charge on 32 3 x 64 lattice. This anal- 
ysis is performed on Fujitsu AP1000 with 1024 
processors. We prepare 830 configurations sepa- 
rated 100 updating sweeps at = 6.00 on 32 3 x 64 
lattice, and then perform blocking twice to obtain 
blocked configurations of size 8 3 x 16. Topological 
quantities are measured using the cooling method 
on 8 3 x 16. Fig. 5 shows the typical cooling history 
of topological charge Q, the sum of the absolute 
value of the local topological density /, and the 
normalized action S, rescaled so that an instan- 
ton configuration has 5 = 1. We use the value 
after 50 cooling sweeps as indicated by the arrow 
in Fig. 5. As for the integrated autocorrelation 
time of the topological charge, r turns out to be 
quite small compared to 100 sweeps which is the 
interval between measurements. In this case we 
have tried to fit p(t) by the exponential function 
of t to get so called exponential autocorrelation 
time, but we can not extract reliable value at this 
stage since p(t) is fluctuating around from t = 
200. 
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